Light-cone coordinates based at a geodesic world line 
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Continuing work initiated in an earlier publication [Phys. Rev. D 69, 084007 (2004)], we construct 
a system of light-cone coordinates based at a geodesic world line of an arbitrary curved spacetime. 
The construction involves (i) an advanced-time or a retarded-time coordinate that labels past or 
future light cones centered on the world line, (ii) a radial coordinate that is an affine parameter 
on the null generators of these light cones, and (iii) angular coordinates that are constant on each 
generator. The spacetime metric is calculated in the light-cone coordinates, and it is expressed as an 
expansion in powers of the radial coordinate in terms of the irreducible components of the Riemann 
tensor evaluated on the world line. The formalism is illustrated in two simple applications, the first 
involving a comoving world line of a spatially-flat cosmology, the other featuring an observer placed 
on the axis of symmetry of Melvin's magnetic universe. 



I. INTRODUCTION 

We continue a research program initiated in Ref. [H, 
which aims to construct and exploit light-cone coordi- 
nates based at an arbitrary world Hne 7 of an arbitrary 
curved spacetime. The retarded coordinates introduced 
in Ref. 0] are denoted {u, r, 9, <j)) and are adapted to the 
future light cone of each point z on the selected world 
line. In this paper we construct advanced coordinates 
{v, r, 9, (j)) which are instead adapted to the past light cone 
of each point z on the world line; for simplicity we take 
the world line to be a geodesic of the spacetime. (The 
advanced coordinates were introduced briefly in Ref. 
this paper provides details that were not given in the 
earlier work.) We collectively denote the light-cone coor- 
dinates by {w, r, 9, <p), with w standing for either u or w 
depending on the context. In both cases the null coordi- 
nate w is constant on each light cone, and it agrees with 
proper time r at the cone's apex. The radial coordinate 
r is an affine parameter on the cone's null generators, 
and it measures the distance away from the world line. 
The angular coordinates 9^ — {9, (f) are constant on each 
one of these generators. The geometrical meaning of the 
light-cone coordinates is clear, and this is one of their 
main virtues. 

The formalism developed in Ref. jl| and in this paper 
incorporates ideas formulated many years ago by Bondi 
and his collaborators 0, Q , and it complements a line of 
research that was initiated by Syngc 5] and pursued by 
Ellis and his collaborators 8, 9, 10, 11] in their work 
on observational cosmology. While the central ideas ex- 
ploited here are the same as with Synge and Ellis, our im- 
plementation is substantially different: While Synge and 
Ellis sought definitions for their optical or observational 
coordinates that apply to large regions of the spacetime, 
our considerations are limited to a small neighborhood 
of the world line. 

The introduction of retarded coordinates was moti- 
vated by the desire to construct solutions to wave equa- 
tions for massless fields that are produced by a pointlike 
source moving on the world line. The retarded coordi- 
nates naturally incorporate the causal relation that exists 



between the source and the field, and for this reason the 
solution takes a simple explicit form (in the neighborhood 
in which the coordinates are defined). The introduction 
of advanced coordinates is motivated instead by the de- 
sire to construct solutions to the Einstein field equations 
that describe black holes placed in a distribution of mat- 
ter or in a tidal environment. Such an application was 
described in Ref. |2] , in which the metric of a tidally dis- 
torted black hole was presented in advanced coordinates. 
In a companion paper we use the guidance offered by 
the advanced coordinates to formulate a light-cone gauge 
for black-hole perturbation theory, and to calculate the 
metric of a black hole immersed in a uniform magnetic 
field. 

A quasi-Cartesian version of the advanced coordinates 
is introduced first in Sec. II B, after reviewing some nec- 
essary geometrical elements in Sec. II A. The metric ten- 
sor in advanced coordinates is constructed gradually in 
Sees. II C through F, and its quasi-Cartesian form is dis- 
played in Eqs. lfT?7|l - lfT^ . 

In Sec. Ill A we combine the results obtained in Sec. II 
with the earlier results of Ref. T] and present the met- 
ric in a general form suitable for both advanced and re- 
tarded coordinates. At this stage the metric is expressed 
in terms of the Riemann tensor evaluated on the world 
line 7. In Sec. Ill B we begin to refine the form of the met- 
ric by decomposing the Riemann tensor into its Weyl and 
Ricci parts, and by further decompositing the Weyl and 
energy-momentum tensors into their irreducible compo- 
nents. This leads us, in Sec. Ill C, to introduce tidal and 
matter potentials that make the basic building blocks of 
the metric tensor. The potentials are displayed in Ta- 
ble I, and the refined form of the metric is displayed in 
Eqs. (pn^ - ITT^ . 

In Sec. IV we carry out a transformation of the met- 
ric from the quasi-Cartesian coordinates x° to the quasi- 
spherical coordinates {r,9^). The transformation, intro- 
duced in Sec. IV A, is the familiar one from flat space- 
time: = rVL°'{9^), or more explicitly, x — r sa\9 coscj), 
y — r siii9 sincj), and z = rcos9. This final expression 
for the metric tensor, in the coordinates {w, r,9^), is dis- 
played in Eqs. (|4.9|) - (|4.12|) . It involves the angular com- 
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ponents of the tidal and matter potentials introduced in 
Sec. III. As shown in Table IV (with the results derived 
in Sees. IV B and C), these are naturally expressed as 
expansions in scalar, vector, and tensor harmonics. The 
required spherical-harmonic functions are listed in Tables 
II and III. 

In Sec. V we present two simple applications of the 
light-cone coordinates. In Sec. V A we apply the formal- 
ism to a comoving world line of a spatially-flat cosmology. 
In Sec. V B we examine the metric near the axis of sym- 
metry of Melvin's magnetic universe [TsL 113. ITsI l . 

Throughout the paper we work in geometrized units 
(G = c = 1) and adhere to the conventions of Misner, 
Thorne, and Wheeler [l^ . 

II. ADVANCED COORDINATES 

The presentation in this Section follows very closely 
Sec. II of Rcf. P]. The material is very similar, but be- 
cause of important differences of sign that occur in var- 
ious places, we present here the details that are specific 
to the advanced coordinates. Other details are omitted 
and can be obtained from Ref. |l| . 



A. Geometrical elements 

We first introduce some geometrical elements on the 
world line 7 at which the advanced coordinates are based. 
The world line is described by parametric relations (t) 
in which t denotes proper time. Its normalized tangent 
vector is = dz^/dr, and we assume that this sat- 
isfies the geodesic equation Du^^/dr — 0. The world 
line is therefore a geodesic of the curved spacetime, and 
this assumption represents a loss of generality relative 
to the construction of retarded coordinates presented in 
Ref. 0] . While it would be a simple matter to restore this 
level of generality, we refrain from doing so in this work. 
Throughout we use Greek indices /i, A, p, etc. to refer 
to tensor fields defined, or evaluated, on the world line. 

We install on 7 an orthonormal tetrad that consists 
of the tangent vector and three spatial vectors e(^. 
These are parallel transported on the world line, so that 
De'^/dT = 0. It is easy to check that this is compati- 
ble with the requirement that the tetrad (u^,e^) be or- 
thonormal everywhere on 7. 

From the tetrad on 7 we define a dual tetrad (e^, e°) 
with the relations e° = — and e° = S°'''g^^el. The 
dual vectors ej^ also are parallel transported on the world 
line. The tetrad and its dual give rise to the completeness 
relations 

(2.1) 

for the metric and its inverse evaluated on the world line. 



The advanced coordinates are constructed with the 
help of a null geodesic segment that links a given point x 
to the world line. This geodesic segment must be unique, 
and we thus restrict x to be within the normal convex 
neighborhood of 7. We denote by (3 the unique, future- 
directed null geodesic segment that goes from x to the 
world line, and x' = z{v) is /3's point of arrival on the 
world line; v is the value of the proper-time parameter at 
this point. To tensors at x we assign the Greek indices 
a, /3, 7, S, etc.; to tensors at x' we assign the indices a', 
/?', 7', S', and so on. 

From the tetrad (m" , e" ) at x' we obtain another 
tetrad (eg,e") at x by parallel transport on /3. By rais- 
ing the frame index and lowering the vectorial index, 
we obtain also a dual tetrad at x: = —QapeQ and 
6a — 3"'^9af3e^- The metric at x can then be expressed 
as 

5a/3 = -e^ie^ + Satelel (2.2) 

and the parallel propagator |^ (also known as the bivec- 
tor of geodetic parallel displacement from x' to x is 
given by 

g"'„,{x,x') = -e^Ua' + e"el,, 

(2.3) 

giix',x) = u^'el + e'^el. 

This is defined such that if A"" is a vector that is par- 
allel transported on /3, then A°'{x) — g°'^,{x,x')A°' [x') 
and {x') — g°'^{x' ,x)A°'{x). Similarly, if pa is a dual 
vector that is parallel transported on /?, then Pa{x) — 
g°'^{x\x)pa'{x') and Pa'{x') = g°'^,{x,x')pa{x). 

The last ingredient we shall need is Synge's world func- 
tion cr (z, x ) la] (also known as the biscalar of geodetic in- 
terval [l3|)- This is defined as half the squared geodesic 
distance between the world- line point z(r) and a neigh- 
boring point X. The derivative of the world function with 
respect to z'^ is denoted iTi_i{z, x); this is a vector at z (and 
a scalar at x) that is known to be tangent to the geodesic 
linking z and x. The derivative of a{z, x) with respect to 
x" is denoted (Tq(z,x); this vector at x (and scalar at z) 
is also tangent to the geodesic. We use a similar notation 
for multiple derivatives; for example, a^a = Vq-V^ct and 
aa(3 = VpVaC, where Vq. denotes a covariant derivative 
at X while indicates covariant differentiation at z. 

The vector — (t^(z, x) can be thought of as a separation 
vector between x and z, pointing from the world line to 
X. When x is close to 7, —a^^{z,x) is small and can be 
used to express bitensors in terms of ordinary tensors at 
z For example, 

cTf^i^ = g,.i^ - ^Rf.xi.pcr^cr'' + ■■■, (2.4) 

fTpa = -g"a(gM- + ^i?MA.pfTV'' + •••), (2.5) 

where g'^^ = g^^{z,x) is the parallel propagator and 
Rp,\iyp is the Riemann tensor evaluated on the world line. 
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B. Definition of the advanced coordinates 

In their quasi-Cartesian version, the retarded coordi- 
nates are defined by 

(2.6) 

The last statement indicates that x' = z{v) and x are 
hnked by the null geodesic segment (3, and we demand 
that this be future-directed from x to x' . 

From the fact that a" is a null vector we obtain 

r:=(^,fax"£'')i/2 = -u„.a"', (2.7) 

and r is a positive quantity because cr" is a future- 
directed vector. We will see in Sec. II C that —r is an 
affine parameter on (3; this property adds credibility to 
the idea that r is meaningful measure of distance from x 
to x' = z{v). 

Another consequence of Eq. 1231 is that 

d"' =r(u"' -f^-'e^'), (2.8) 
where $7° := a;°/r is a frame vector that satisfies 

A straightforward calculation reveals that under a dis- 
placement of the point x (which induces a displacement 
of x'), the advanced coordinates change according to 

dv = -ladx", (2.9) 
dx" = -el,a°''p,u^' dv - el,a°''i:sdx'^ , (2.10) 

where la '■= — Ca/r is a future-directed null vector at x 
that is tangent to (3. 

C. Advanced distance; null vector field 

If we keep x' linked to x by the relation a{x^x') — 0, 
then 

r{x) = -aa'{x,x)u'^\x') (2.11) 

can be viewed as an ordinary scalar field defined in a 
neighborhood of 7. We can compute the gradient of r by 
finding how r changes under a displacement of x (which 
induces a displacement of x'). The result is 

V^jr = [aa'fi'u'^'u^'Ysi ~ ct^'/jw"'. (2.12) 
Similarly, we can view 

= -^^^ (2.13) 

as an ordinary vector field, which is tangent to the con- 
gruence of null geodesies that converge to x' . It is easy 
to check that Eqs. ifTT^ and imply 

rV„r = -l. (2.14) 



In addition, combining the general statement a" = 
~9°'a'<^"' with Eq. gives 

r =5"„,(^"'-r!'^e^'); (2.15) 

the vector at x is therefore obtained by parallel transport 
of w"' - f]°e^' on (3. From this and Eq. we get the 
alternative expression 

r = e^-n''e^, (2.16) 

which confirms that is a future-directed null vector 
field (recall that fl'^ = x"" /r is a unit frame vector). 

The covariant derivative of can be computed by 
finding how the vector changes under a displacement of 
X. This calculation reveals that satisfies the geodesic 
equation in afhne-parameter form, l^Vpl" — 0, and 
Eq. (|2.14() informs us that the affine parameter is in fact 
—r. A displacement along a member of the congruence is 
therefore described by dx" = —I" dr. Specializing to the 
advanced coordinates, and using Eqs. H2.9|l . (|2.1(J|) . and 
(|2.13() . we find that this statement becomes dv = and 
dx"' — (x"/r) dr, which integrate to w = constant and 
— rff^, respectively, with fi" representing a constant 
unit vector. We have found that the congruence of null 
geodesies that converge to x' is described by 

w = constant, = rr2''(6'^) (2.17) 

in the advanced coordinates. Here, the two angles 9^ 
{A — 1,2) serve to parameterize the unit vector fi", 
which is independent of r. 

Finally, we state without proof that I" is hypersurface 
orthogonal (the proof is contained in Ref. (ij). This, to- 
gether with the property that I" satisfies the geodesic 
equation in affine-parameter form, implies that there ex- 
ists a scalar field v{x) such that 

= -Vc-y. (2.18) 

This scalar field was already identified in Eq. 1)2. 9|) : it is 
numerically equal to the proper-time parameter of the 
world line at x' . We conclude that the geodesies to 
which P is tangent are the generators of the light cone 
V = constant. As Eq. H2.17|l indicates, a specific gener- 
ator is selected by choosing a direction il" (which can 
be parameterized by two angles 0^), and —r is an affine 
parameter on each generator. The geometrical meaning 
of the advanced coordinates is now completely clear; the 
construction is illustrated in Fig. 1. 

D. Frame components of tensors on the world line 

The metric at x in the advanced coordinates will be 
expressed in terms of frame components of tensors eval- 
uated on 7. We shall need, in particular, 

RaQhoiv) Ra'Yti'S' V? It* , 

Rabcoiv) ■■= i?a'7'/3'5'e^'efe2'u*', (2.19) 
Racbdiv) ■■= Ra'YP'S' ea'eteb ■ 
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FIG. 1: Advanced coordinates of a point x relative to a world 
line 7. The advanced time v selects a particular light cone, the 
unit vector il" := x'^/r selects a particular generator of this 
light cone, and the advanced distance r selects a particular 
point on this generator. 



These are the frame components of the Riemann tensor 
evaluated on 7; these quantities depend on v only (recall 
that V is numerically equal to the proper-time parameter 
on the world line) . We next form the useful combinations 



Pa, 



- RaObO — RacbO^'^ — RbcaO^^'^ + Racbd^^'^^'^ 
= Pba. (2.20) 
Pa ■= Pab^^ = RaObQ^^ — RabcO^^^'^ , (2-21) 

P := Pa^l" ^ Ramo^"^\ (2.22) 

in which the quantities fi" :— x"" /r depend on the angles 
9^ only — they are independent of v and r. 



F. Metric near 7 

It is straightforward to invert the relations of 
Eqs. ^(TT^ and ifT^ and solve for da;" and dx". 
The results are 



l + ]^r^P + 0{r^) 
1 



dv 



- -r\Pa - PVLa) + 0{r^) 


1 



^a 

dx". 



(2.25) 



dv 
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r^{P%~ P''nb) + 0{r^) dx\ (2.26) 



These relations, when specialized to the advanced coordi- 
nates, give us the components of the dual tetrad (e° , e°) 
at x. The metric is then computed by involving the com- 
pleteness relations of Eq. (|2.1|l . We find 



9vv 
9va 

gab 



= -1-r^P + 0{r^), 

^r^{Pa-Pna)+0{r^), 



Sab 



-yiPab 



Pn^h 



-naPb + P^a^b)+0{r''). 



(2.27) 
(2.28) 

(2.29) 



We see that the metric possesses a directional ambiguity 
on the world line: The metric at r = still depends on 
the vector Vt"^ := x°'/r that specifies the direction to the 
point X. The advanced coordinates are therefore singu- 
lar on the world line, and tensor components cannot be 
defined on 7. This poses no particular difficulty because 
we can always work, as we have been doing, with frame 
components of tensors instead of tensorial components. 



III. LIGHT-CONE COORDINATES; 
DECOMPOSITION OF THE RIEMANN TENSOR 

A. Retarded and advanced coordinates 



E. Coordinate displacements near 7 

The changes in the quasi-Cartesian advanced coordi- 
nates under a displacement of x are given by Eqs. (|2.9() 
and (|2.1()|l . In these we substitute the expansions for 
CTq,'/3' and aa'p that appear in Eqs. (|2.4|1 and l|2.5|l . as 
well as Eqs. (|2.8f) and (|2.1t)|) . After a straightforward 
calculation, we obtain the following expressions for the 
coordinate displacements: 



dv — 
dx" = 



(e°dx") 
2 

1 



+ fla(ei dx-^) 



(2.23) 



0(r3) 



5\ 



-r\P\ + 2P-'n,) + 0(r3)J (e^ d:r").(2.24) 

Notice that the result for dv is exact, but that dx'^ is 
expressed as an expansion in powers of r. 



The developments of Sec. II parallel very closely the 
construction of retarded coordinates described in Ref. j3] ■ 
The combined set of results is a coordinate system {w, x°') 
that refers either to past light cones (advanced coordi- 
nates, w = V, 7] = +1) or to future light cones (retarded 
coordinates, w = u, rj = —1) centered on a geodesic world 
line. In either case the metric is expressed as 



with 



9ww 
V9wa 

gab 



Pa 



-1 
^a 

Sab 



r^P + 0{r^), 
2 
3 

fla^b 



pn„ 



1 



Pn 



-Pa^b 



-naPb + pnaiib) + 0{r''), 

RaObO ~ VRacbO^"^ ~ VRbcaO^'^ 



(3.1) 
(3.2) 

(3.3) 
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Pa 
P 



+ Racbd^'^^ 
RaObO^^ — r]RabcO^^^'^ 
RaObQ^'^^^ ■ 



(3.4) 
(3.5) 
(3.6) 



The Riemann tensor is evaluated at the ad- 
vanced/retarded point x' = z{w), and its frame 
components are defined as in Eqs. (|2.19l) . Our subse- 
quent developments will apply to the general metric 
of Eqs. (|3.1II - (|3.3|I . They will not distinguish between 
the advanced and retarded coordinates, and we will 
collectively refer to them as light- cone coordinates. 
As was stated previously, w stands for either v or m, 
and T] is an indicator that takes the value +1 for the 
advanced coordinates, and the value —1 for the retarded 
coordinates. 



B. Decomposition of the Riemann tensor 

To bring the metric to a more explicit form we decom- 
pose the Riemann tensor into its Weyl and Ricci parts, 
and we involve the Einstein field equations to relate the 
Ricci tensor to the energy-momentum tensor of the mat- 
ter distribution. This gives 

■^9a'[Y9s']f3' T, (3.7) 

where Ca'iS'-y'S' is the Weyl tensor, Ta'p' is the energy- 
momentum tensor, "^T := T'^, is its four-dimensional 
trace, and where the square brackets indicate antisym- 
metrization of the enclosed indices. 

We next project the Weyl tensor onto the tetrad 
(m" , e" ) and decompose the projections into irreducible 
components, according to 



CaObO 


— ^ab, 


CabcO 


— ^abpB^c 


Gabcd 


= Sac£bd 



(3.8) 



<^bc£ad — ^ad^bc + <^bd£ac, 



where Sabc is the flat-space permutation symbol. The 
electric components of the Weyl tensor are denoted £ab, 
while the magnetic components are denoted Bab- These 
tensors are symmetric and tracefree, so that, for example, 
Sba = £ab and B"^^ — 0. Because the Weyl tensor is 
evaluated on the world line, £ab and Bab are functions of 
the null coordinate w only. 

We perform similar operations on the energy- 
momentum tensor, and introduce the notation 



Too 

Toa 
Tab 



Sab - 



(3.9) 



-5abT, 



where Sab is symmetric and tracefree, and T :— S Tab = 
T„'/3'(5"^e^'ef ) = T„,0,(g"'^' + u"'uf^') = '^T + p is the 



TABLE I: Tidal and matter potentials. Each potential is iden- 
tified with a sans-serif superscript that specifies its multipole 
content. A potential labeled with a "q" is a quadrupole field, 
and one labeled with an "d" is a dipole field. The vector and 
tensor potentials are all orthogonal to fi" , and all tensors are 
symmetric and tracefree. 

r = Scd^'^n'' 

"Sab = 2(5,'= - f7an=)(5(,'' - nb^'')£cd + {dab - fiafib) 

'Bab = Eap.nPB'^.iS^ - n'^Qb) + ebp,nPB'',{S% - Q'^Qa) 
% = {6,^ - n„n=)j, 

SS = Scd^^^^ 

SSa — {Sa — ^la^'')Scd^'^ 



spatial trace of the energy- momentum tensor. The quan- 
tity p represents the mass-energy density measured by 
an observer moving on the world line 7, ja is the flux 
of mass-energy traveling in the direction of the base vec- 
tor e" , Sab is the tracefree part of the stress tensor, and 
ir is an isotropic pressure. These quantities also are 
functions of w only. 

Substituting Eqs. (pD^ and into Eq. iPTTjl . and 

then this into Eq. (|3.4|) , produces a decomposition of Pab 
into its irreducible pieces. We obtain 

Pab — 2£ab — ^a£bc^'^ — ^b£ac^'^ + 5ab£cd^'^^'^ 

-r,{eapgnPBP, + ebpgnPBP^) 
+ y(35ab-2f^af^b)p 

- ATTTjljaflb + jb^a - ^Sabjc^") 

- 47r(f2aS'bcf^'^ + ^bSac^"^ — SabScd^'^^'^) 

+ Y^abT. (3.10) 



C. Tidal and matter potentials 

At this stage it is useful to involve the irreducible quan- 
tities £ab{w) and Bab{w) in the definition of a number of 
tidal potentials. We also involve ja(w) and Sabiw) in the 
definition of matter potentials. These potentials, which 
are displayed in Table 1, form the elementary building 
blocks of the metric tensor. Each potential is identified 
by its multipole content. For example, ''£ := £ab^'^^^ 
is a quadrupolar potential by virtue of the fact that 
£ab is symmetric and tracefree. As another example, 
"1? ■— ja^^" is a dipolar potential. Table I also intro- 
duces vectorial and tensorial potentials that possess the 
property of being transverse, meaning that each vector or 
tensor is orthogonal to the unit frame vector fl'^. Finally, 
the tensor potentials ''£ab and ''Bab have the additional 
property of being symmetric and tracefree. 



6 



It is easy to check that Pab is expressed in terms of the 
tidal and matter potentials as 

Pab = 'Sab + 2n(^a'Sb) + na^b'S 
-rj{'Bab + 2n^a^Bb)) 

+ 4Tr[Sab - i^a^b)p + -Y^a^bP 

+ STr7][{Sab - ^ab)'] - ^ia'jb)] 

+ 4T:{6ab - ^a^b)SS - STTflf^aSSb) " ATrVlanbSS 

(3.11) 



— {Sab - fla^b)T + —ila^bT. 



We observe that Pab is now decomposed into transverse- 
transverse components that are fully orthogonal to Qa, 
transverse-longitudinal components that are partly or- 
thogonal to and partly aligned with f2a, and longitudinal- 
longitudinal components that are proportional to ^a^b- 
Contracting Eq. (|3.11() with fl^ produces 



IV. ANGULAR COORDINATES 
A. Transformation to angular coordinates 

Because the frame vector fi" := x°'/r satisfies 
Sab^"^^^ = 1, it can be parameterized by two angles 9^. 
A canonical choice for the parameterization is 



$7" = (sin 6 cos (j>, sin 6 sin 0, cos I 



(4.1) 



It is then convenient to perform a coordinate transforma- 
tion from to {r,9^) using the relations x° = ril,°'{9^). 
(Recall from Sec. II C that the angles 9^ are constant 
on the generators of the light cones w = constant, and 
that ±r is an affine parameter on these generators. The 
relations = rfi" therefore describe the behavior of 
the generators.) The differential form of the coordinate 
transformation is 



rn'X d9^ 



(4.2) 



Pa = '^Ea + ^a'S-rBa + Y^ap-^^v'ja 

-A^SSa-A^^aSS+^naT, (3.12) 



and contracting this with gives 



47r 47r 
P=''£ + —P~A'kSS+—T. (3.13) 



Substituting Eqs. l|TTT|l - l|ni|l into Eqs. 
produces our final expression for the metric tensor in the 
quasi-Cartesian version of the light-cone coordinates. We 
obtain, after simplification. 



47r 

9n,n, = -l-r2r-yr2(p-355 + r)+0(r3),(3.14) 



where the transformation matrix 

satisfies the identity fla^A = 0- 
We introduce the quantities 

^AB ■= Sab^A^B^ 



(4.3) 



(4.4) 



which act as a (nonphysical) metric on the submanifold 
spanned by the angular coordinates. In the canonical 
parameterization of Eq. H4.1|) . Uab = diag(l, sin^ 0), and 
the metric is that of a round two-sphere of unit radius. 
We use the inverse of ^ab , denoted , to raise upper- 
case Latin indices. We then define the new object 



(4.5) 



which satisfies the identities 



+ ^r^r'{SSa + V%)+0{r'), 

gab = 5ab ~ ^a^b - ^r^i'Sab ~ V'^^ab) 



An 



3 

+ 0{r^) 



(3.15) 



{Sab - ^a^b) (p + 2ri "j + SS+ |T) 



(3.16) 



^a^B = ^B^ 



^A^b 



s^ — n'^rth 



(4.6) 



The first result is a direct consequence of the definition, 
and the second result follows from the fact that both sides 
are symmetric in a and b, orthogonal to fla and fi^, and 
have the same trace. 

The Levi-Civita tensor on 5*^ is constructed as 



£AB ■— Sabc^A^^B^'^' 



(4.7) 



We observe that the metric is neatly expressed in terms 
of the monopolar "potentials" p and T, the dipolar po- 
tentials ''j and ''ja, the quadrupolar tidal potentials 
''Sa, ''Sab, ''Ba, ''Bab, and the quadrupolar matter poten- 
tials SS and SSa- We observe also that g^a contains 
both longitudinal and transverse pieces, while gab is fully 
transverse. We recall that the coordinates are advanced 
when r/ = +1 (then w = v) and that they are retarded 
when r/ = — 1 (then w = u). 



where Sabc is the Cartesian permutation symbol; in the 
canonical coordinates we have e^^ = sin 9. 

We let Da denote the covariant derivative operator 
compatible with ^ab, so that Da^bc =^ 0. It is easy to 
show that DaSbc = and 



Ms- 



(4.8) 



When we apply the coordinate transformation of 
Eq. to the metric of Eqs. we find that 
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the only nonvanishing components of the metric tensor 
are now given by 

Att 

g^un, = -l-r^''£-—r^{p-3SS + T)+0{r^),i4.9) 

(4.10) 



IwA = - 



^T^r^SSA + v'3A)+0{r^), (4.11) 

QAB = r^^AB - ^^r^i^SAB - V^Bab) 

- —r^{lAB{p + 2r,'j + SS + \T) 



+ 0(r5), 



(4.12) 



where "Ea ■= "Sai^A^ "Sab ■= "Sabi^A^'h, and so on. 
The results g^jr = drr = 0, and grA = are exact, 
and they follow from the light-cone nature of the coor- 
dinates. [For example, for the advanced coordinates we 
have Ic, = (-1,0,0,0) and 1°' = (0,-1,0,0), where we 
use the ordering {v,r,9,(f>); these relations imply that 
gyr — 1 and grr = grO — gr4> = 0-] Once more we recall 
that the coordinates are advanced when 77 = +1 (then 
w = v), and that they are retarded when 77 = —1 (then 
w = u). 



B. Tidal and matter potentials in spherical 
coordinates 

According to Table I, the tidal potential ''£ is defined 
by = Scd^'^^'^- Differentiating this with respect to 6^ 
gives Da ''£ = 2n\£cd^'^- In view of the identity Slcf^^i = 
we may write this as Da''£ = '^^A.i^a - ^a^'')£cdi^'^- 
Referring once more to Table I, we see that this is 
Da "*£ = 2il^ ''£a and we conclude that 



''£a = \da''£. 



(4.13) 



Acting on '^£ with two derivative operators gives 
DaDb''£ = 2£cd^\n% + 2£cd^''n%. Using Eq. 
produces DaDb''£ = 2£cd^\n% - 2£cd^''n'^nAB = 
2£crff^A^^i-2fiAB "£■ We write this in the form {DaDb+ 
SriAs)''^ = 2£,d^\n% + nAB''£ - 2£,rfr!^r!^((S,^ - 
r!af^")('5fc'' - ^b^'^)£cd + ^A^ei^ab - na^b)''£, after in- 
volving Eq. (|4.4|l . Consulting Table I once more, we see 
that the right-hand side is equal to f^^fi^ '^£ab "^ab 
and we conclude that 



"^AB = {DaDb + i^AB) '£■ 



(4.14) 



We observe that this tensor is tracefree, because 
f]^^ "Sab = ri^^r!:^^!^ '£ab = {5"''' - f^"f^') '£ab = 0, af- 
ter involving Eq. and the fact that "fab is transverse 
and tracefree (in the Cartesian sense). The equation 
{n'^^DADB + Q) "£ = 0, which we obtain from Eq. g33||, 
reveals that ''£{w, 9^) is a spherical-harmonic function of 
degree 1 = 2. 



We define a magnetic potential "i? := Bcd^'^^'^ and 
differentiate it with respect to 9^, giving Db'^B = 
2ri^;Bcdf^''- We next multiply this by the Levi- 
Civita tensor of Eq. l|lT|l and get -EaDb^B = 
~2eapq^A^P^^''^B^cd^'^ ■ Using Eq. gS]) and the an- 
tisymmetry property of the permutation symbol, this is 
-e/Db^B = -2n'Xeapqfl'^B''^n'' = 2VlYBa, and we 
conclude that 



^6. 



1 



e/Db^B, "B-.^Bcd^^n''. 



(4.15) 



Acting on "S with two derivative operators and multi- 



plying by the Levi-Civita tensor gives —e^ DbDc^B = 
2n'^eapq^^B'^f^n'^g + 2eAB "S- Symmetrizing with respect 
to A and B and consulting Table I yields —{s^Db + 
e§DA)Dc "B = 2f}^fi|j '<Bab, and we conclude that 



'Bab = --{e^DB+e§DA)Dc'B. 



(4.16) 



This tensor is tracefree, because 51 "Sab = 
—e^'-^ DbDc'^B = by virtue of the symmetry of 
DbDc'^B and the antisymmetry of the Levi-Civita ten- 
sor. 

Similar results can be obtained for the matter po- 
tentials. Differentiating "'j := ja^'^ produces Da'^J — 

ja^'A — jc{5a ~ ^a^'^)^'A = ^A'l?a, and WC COUCludc 

that 



Finally, 



'jA = DA']. 



SSa — -DaSS 



(4.17) 



(4.18) 



follows after a calculation similar to the one leading to 
Eq. mUl . 



C. Decomposition in spherical harmonics 

The results obtained in the preceding subsection in- 
dicate that the tidal and matter potentials that appear 
in the metric of Eqs. H4.9|l - H4.12|l can all be obtained by 
covariant differentiation of the scalar potentials "5, "B, 
■y, and SS. These are functions of the null coordinate 
w and the dependence on the angles 9^ appears in the 
factors n''{9^). 

This angular dependence can be made more explicit 
by involving spherical-harmonic functions. Let 



|yl,0 yl.lc 



(4.19) 



be a set of real, unnormalized, spherical-harmonic func- 
tions of degree 1 = 1. And let 

y2m _ |y2,0 y2,lc y2,ls y2,2c y2,2s| (4 20) 

be a set of real, unnormalized, spherical-harmonic func- 
tions of degree 1 = 2. The abstract index m describes 
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TABLE II: Scalar and vectorial harmonics of degree I — 1, labeled by the abstract index m = {0, Ic, Is}. The odd-parity 
vectorial harmonics X}[" are not required, and the tensorial harmonics Y}]^ and X]{^ vanish identically. 



m 









Ic 


Is 


y Im 




cos 9 




sin 9 cos 


sin 9 sin 


T^lm 




— sin 9 




cos S cos <j) 


cos 9 sin 


TAlm 









— sin ^ sin 


sin cos 


TABLE III: Scalar, vectorial, and tensorial harmonics of degree / — 


2, labeled by the abstract index 


m = {0, Ic, Is, 2c, 2s}. 


m 





Ic 


Is 




ZS 




-(Scos^e- 1) 


2 sin 9 cos 9 cos (j) 


2 sin 9 cos ^ sin 


sin^ 9 cos 20 


sin^ 9 sin 20 




6 sin 8 cos 6 


2(2cos^6'- l)cos(?!) 


2(2cos^6l - l)sin(?; 


> 2 sin S cos S cos 20 


2 sin 9 cos ^ sin 20 







—2 sin ^ cos 9 sin 


2 sin 6* cos 9 cos 


-2sin^6'sin20 


z sm t7 cos Z(p 


^ee 


-Ssin^e 


—2 sin 9 cos S cos 


—2 sin 9 cos ^ sin 


(cos^ e + 1) cos 20 


(^cos t7 + 1 j sm z(p 


V2m 





2 sin^ 6^ sin (f> 


— 2 sin^ ^ cos (f> 


— 2 sin 6* cos 9 sin 20 


2 sin 6* cos S cos 20 


V2m 


3 sin** 9 


2 sin'^ 6* cos 9 cos 


2 sin^ 61 cos 9 sin 


-sin^ 6l(cos^ 6* + 1) cos20 - 


sin^ ^(cos^ e + 1) sin 20 


\>'2m 

^0 





2 cos 9 sin 


—2 cos 9 cos 


2 sin ^ sin 20 


—2 sin 9 cos 20 


\^2m 
-^4, 


6 sin^ 9 cos 9 


2sin6'(2cos2 6l - l)cos<;/) 


2 sin 9(2 cos^ S - 1) sin 2 sin^ 9 cos S cos 2<;!) 


2 sin^ 9 cos 6* sin 20 


-^00 





—2 sin sin 


2 sin cos 


2 cos sin 20 


—2 cos cos 20 


V"2m 
-^04, 


-Ssin^e 


—2 sin^ S cos ^ cos (f> 


—2 sin^ ^ cos S sin (p 


' sin ^(cos^ ^ + 1) cos 20 


sin0(cos^ + 1) sin 20 


"V2m 

^4>4 





2 sin^ 61 sin 


-2sin^ 6'cos(/) 


-2 sin^ 6* cos 61 sin 20 


2 sin^ 9 cos 61 cos 20 



Apart from notation and normalization, these definitions 
agree with those of Regge and Wheeler [T9| . We note that 
the tensorial harmonics and are symmetric and 
tracefree. 

The decompositions of the vectorial and tensorial po- 
tentials in terms of vectorial and tensorial harmonics are 
displayed in Table IV. They are obtained by substituting 
Eqs. (jOT|| - ||i:^ into Eqs. ^T^-^^. 

The most explicit form for the metric tensor is obtained 
after substituting the spherical-harmonic decompositions 
of Table IV, along with the spherical-harmonic functions 
listed in Tables II and III, into Eqs. This 
leads to long expressions, but in practical applications 
it may happen that only a few frame components among 
S'ob^ Sab, P, ja, Sab, and T are nonzero; in such cases only 
a few harmonic components among ''Sm, ''Bm, "'jm, and 
SSm will contribute to the metric, and the expressions 
will simplify. We shall encounter such cases in the next 
section. 



V. APPLICATIONS 

A. Comoving observer in a spatially-flat cosmology 

To illustrate how the formalism works we first consider 
the world line of a comoving observer in a cosmological 
spacetime. The global metric is 

ds2 ^ -dt^ + a^{t){dx^ + dy^ + dz^) , (5.1) 



the dependence of the spherical-harmonics on the angle 
0; the numerical part of the label refers to the azimuthal 
index m, and the letter indicates whether the function is 
proportional to cos(m0) or sin(rn0). Explicit expressions 
are listed in Tables II and III. 

We decompose the scalar potentials according to 





6^) 


m 


e% 


(4.21) 






m 




(4.22) 






m 




(4.23) 


SS{w, 






{0% 


(4.24) 



m 



in terms of their harmonic components '*£m, '*Bm, '^jm, 
and SSm ■ These are in a one-to-one correspondence with 
the frame tensors £ab{w), Bab{w), ja{w), and Sabiw); the 
relationships are displayed in Table IV. 

The derivatives of the scalar potentials will be ex- 
pressed in terms of derivatives of the spherical-harmonic 
functions. The vectorial harmonics are 

rj'^ = DaY^"", = -efDBY^"", (4.25) 

and the tensorial harmonics are 

r'^nAB, Y'^B^[DADB + \l{l+mAB\Y"^ (4.26) 

and 

X% = -\{e^DB + e§DA)DcY''^. (4.27) 
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TABLE IV: Spherical-harmonic decomposition of the tidal and matter potentials. In the first part of the table we list the 
definitions of ''jm{w), ''£^{'w), ''Bm{w), and SSr,^{w) in terms of ja{w), Eabiw), Bab{w), and Sab{w). In the second part of the 
table we display the spherical-harmonic decompositions of the tidal and matter potentials. 



'jo = ja 

'jlc = jl 
'jls = j2 


ro = i(£ii + £:22) 

"^Eis ~ £23 

''£2c = ^(£^11 ~ £22) 
''£2s = £12 


-Bo = |(Bii +S22) 

"Bic = Bi3 

"Bi^ = S23 

''B2C = i(Bii -S22) 

■'B2« = S12 


SSo ~ ^(•S'li + S22) 
SSic ~ S13 

SSls = ^23 

SS2C ~ ^('S'li — 522) 

SS2S ~ Si2 


'j = %Y''^ 

d • d ■ \/-lm 
J A = Z^n, J'^'^A 


= Er. "i^my 

^fA = iE.^£-y|"^ 

'Sab = E. "fmFiS 


'Ba = i E. 

'Bab = E. ^Sr.Xl'S 


ss^j:^ ss^y""^ 

SSa = ^ Em SSmY^"' 



where a{t) is an arbitrary scale factor; for simplicity we 
take the cosmology to be spatially flat. This application 
was already presented in Ref. jj, but we generalize it 
here from the retarded coordinates considered there to 
Hght-cone coordinates of both types (retarded and ad- 
vanced). Fm'thermore, the decomposition of the energy- 
momentum tensor into irreducible parts was not accom- 
plished in the earlier paper, and this decomposition adds 
insight to our earlier results. 

Without loss of generality we take our observer to be 
at the spatial origin of the global coordinate system {x — 
y = z — 0), and her velocity vector is given by 

= (1,0,0,0) (5.2) 

in the ordering (t, x, y, z) of the cosmological coordinates. 
This vector satisfies the geodesic equation, and t is proper 
time for the observer. We wish to transform the metric 
of Eq. (|5.1II to light-cone coordinates {w,r,9^) centered 
on the world line of this observer. 

To do so we must first construct a triad of orthonormal 
spatial vectors e^. A simple choice is 

= (0,a-\0,0), 

= (0,0,a-\0), (5.3) 
= (0,0,0,a-i); 

these vectors are all parallel transported on 7. 

According to Eq. H3.8I) and the fact that the Weyl ten- 
sor of the spacetime vanishes, we have 

Sab = Bab = 0. (5.4) 

And according to Eq. (|3.9() and a simple computation, 
we have ja = Sab — and 

p=|-(d/a)^ T^~^[2-d/a+ia/af]. (5.5) 

Here the scale factor is expressed in terms oi w = 
[proper time on 7] by simply making the functional sub- 
stitution a(t) — s- a(w); overdots indicate differentiation 
with respect to w. Recall that p is the mass-energy den- 
sity measured by the observer, and that is the mea- 
sured pressure of the cosmological fluid. 



The vanishing of £ab, Bab, ja, and Sab implies that the 
metric is spherically symmetric around 7 (this does not 
come as a surprise). After substituting Eqs. H5.4I) and 
(|5.5|) into Eqs. H4.9|) - (|4.12|l . a short calculation reveals 
that the metric components are given by 

grow = -l + r^{ii/a) + 0{r^), (5.6) 
9wr = V, (5-7) 
9wA = O(r^), (5.8) 

gAB = r^nAB^l + ^r^[a/a-{k/af]+0{r^)y{5.9) 

We recall that the scale factor and its derivatives are 
functions of the null coordinate w. When the scale factor 
behaves as a power law, a(t) oc with a a constant, we 
have d/a = —a{l — a)/w'^ and d/a — {d/aY = —a/uP'. 
When instead the scale factor behaves as an exponential, 
a{t) oc e^* with H a constant, we have d/a — and 
d/a — {d/aY — 0- 

B. Static observer in Melvin's magnetic universe 

Melvin's magnetic universe 0, [T3 . ITU is a static, 
cylindrically symmetric spacetime that is filled with a 
magnetic field held together by gravity. The exact so- 
lution to the Einstein-Maxwell equations that describes 
this situation consists of the metric 

and the vector potential 

A" = ^BAifi", (5.11) 

where ip" = dx" /dip is the spacetime's azimuthal Killing 
vector. We have introduced 

A:=l + iBV, (5.12) 

and the constant B measures the strength of the mag- 
netic field. The metric and the vector potential are ex- 
pressed in cylindrical coordinates {t, p, z, (p). 
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The metric of Eq. (|5.10() can be decomposed in terms of 
a tetrad of orthonormal vectors. We introduce a "Carte- 
sian" frame described by 



(A-\0,0,0), 

(O, A^"^ cos iy9, 0, —Ap~^ sin i^) , 
(O, A^^ sin If, 0, Ap^^ cos (p) , 
(0,0,A-i,0). 



(5.13) 
(5.14) 
(5.15) 
(5.16) 



It is easy to check that the inverse metric can be ex- 
pressed as g 



af3 _ 



-0 ^0 



„/3 



13 



2 



feg. It is also 



easy to check that in this tetrad, the electromagnetic field 
tensor has 



-B3 :— F12 :— i^Q/3e"e2 — -r^ 



B 
A2 



(5.17) 



as its only nonvanishing component. 

We wish to consider a static observer in Melvin's mag- 
netic universe. To ensure that this observer moves on a 
world line 7 that is a geodesic of the spacetime, we place 
him on the axis of symmetry at p — Q. The observer has 
a velocity vector given by = Cg (p = 0) , and e" (p = 0) 
is a triad of parallel-transported vectors on 7. 

A straightforward computation reveals that the metric 
of Eq. H5.10|l comes with a Weyl tensor whose nonvanish- 
ing electric components are 



^11 — £2 



1 



-B' 



33 



-B' 



(5.18) 



The magnetic part of the Weyl tensor vanishes: Bab — 0. 
A computation of the energy-momentum tensor (either 
from the metric or from the electromagnetic field tensor) 
reveals that 



(5.19) 



Sii — S2 



B^ 
127r' 



S33 — 



B^ 



and 



B^ 

Stt' 



(5.20) 



(5.21) 



while ja = 0. (RecaU that T := S'^^Tab is the three- 
dimensional trace of the energy-momentum tensor; the 
four-dimensional trace is = T — p, and it vanishes 
by virtue of the conformal invariance of Maxwell's equa- 
tions.) These relations imply that ''Eq — ^B^ and 
SSq — i?^/(127r) are the only nonvanishing harmonic 
components of the tidal and matter potentials. 

Making the substitutions from Eqs. H5.18|l - H5.21(l . Ta- 
ble III, and Table IV into Eqs. ^^-^^J^, we find that 
the metric components in light-cone coordinates are 



Qww 

9w0 



I - -B'^r'^ sin^ 
2 



1 



rjB-'r-' sm6cos6 + Oir^), 



6 



B^r'^ sin^ 6* + 0(r5). 



(5.22) 
(5.23) 

(5.24) 
(5.25) 



2 • 2 
r sm ( 



6 



B^r'^ si-a^ e + 0{r^). (5.26) 



As expected, the metric is axially symmetric, but the full 
cylindrical symmetry of the spacetime is not revealed by 
the light-cone coordinates. 
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